Universal temperature dependence of the conductivity of a strongly disordered 

granular metal 



A. R. Akhmerov, A. S. Ioselevich*, 
L. D. Landau Institute for Theoretical Physics, 119334 Moscow, Kosygina 2, Russia 

(Dated: February 6, 2008) 

A disordered array of metal grains with large and random intergrain conductances is studied 
within the one-loop accuracy renormalization group approach. While at low level of disorder the 
dependence of conductivity on log T is nonuniversal (it depends on details of the array's geometry), 
for strong disorder this dependence is described by a universal nonlinear function, which depends 
only on the array's dimensionality. In two dimensions this function is found numerically. The 
dimensional crossover in granular films is discussed. 



Introduction. Unusual temperature dependence of 
conductivity a of granular metals has been in the focus 
of both experimental (see, ej^, yl 3, Ij and references 
therein) and theoretical 0, WB, H S 11 EJ H 
studies in recent years. In particular, the logarithmic 
T-dependence of resistivity 

p = A(l - odnT), (1) 

was experimentally established in Q for low-resistivity 
samples, while high resistivity samples, studied by 
the same group, demonstrated the characteristic Efros- 
Shklovskii law p cx exp{(T ss /T) 1 / 2 } (see [HI El for the- 
oretical interpretation). The empiric dependence Q was 
found to be valid in the wide T-range 30-300 K, where 
the resistivity p varied by a factor of order 3. 

A similar weak T-dependence of low-resistivitygran- 
ular films was observed also in Although in jj this 
dependence was interpreted as a power-law p = AT~ a 
with small a — 0.117, it could with the same accuracy 
be fitted by the law (JTJ. 

The dependence can not be explained by the usual 
weak-localization effects because of many reasons. In 
particular, the dependence was apparently found in 
both two- and three-dimensional samples, and, secondly, 
it appears to be insensitive to the magnetic field. 

The logarithmic T-dependence of the conductivity (not 
the resistivity!) was theoretically derived by Efetov and 
Tschersich 4]. It arises due to corrections to conduc- 
tivity, originating from local quantum fluctuations of in- 
tergrain voltages, a la Ambegaokar, Eckern, and Schon 
|14j . Under the conditions go 3> 1 (low resistivity at 
high temperature) and god «T« Ec (moderately low 
temperature), Efetov and Tschersich have found that 
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for an ideal lattice of identical grains with identical inter- 
grain conductances. Here go is the dimensionless (mea- 
sured in the units of e 2 /2irH) intergrain conductance at 



high temperature T ~ Tr = goEc, 8 is the intragrain 
level spacing, z is the coordination number of the lat- 
tice, and Ec 3> 8 is the charging energy of an individual 
grain. Note, that tr = H/Tr is nothing else, but the 
characteristic local charge relaxation time in the system. 
The formula J2J is valid for a lattice of arbitrary dimen- 
sionality. 

Later, in the paper ||| the theory of Efetov and Tsch- 
ersich was generalized for the case of nonideal - random 
- array of metal grains. It was shown, that the system 
of one- loop RG equations, describing the renormalization 
of individual intergrain conductances gij upon lowering 
of temperature, can be written in the form 



din, gij 
dt 



-2R l3 ({g kl }), 
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where t = In (Tr/T) - is the "RG time", and Rij are di- 
mensionless resistances of the entire network between the 
grains i and j. These resistances are some complicated 
nonlocal functions of all the conductances in the array; 
to find them one first would have to solve the set of the 
Kirchhof equations 
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for voltages V k , and then find Rij({gki}) = V i 

(ij) 

Vj ' . This task is analytically not feasible for general 
random inhomogeneous network of conductances. 

For the case of weak disorder it was shown in || that 
the theoretical dependence of a on In T considerably de- 
viates from the linear law only in the vicinity of the 
metal-insulator transition, where the role of disorder be- 
comes crucial. In the major part of the domain of the 
applicability of the one-loop RG approach the Efetov- 
Tschersich law (JSJ remains a good approximation even 
for a moderately strong disorder, when the relative fluc- 
tuations of intergrain conductances are of order of unity. 
At least this is true for the case of square lattice, which 
was studied numerically in ||. 

In the present Letter we are considering the case of 
extremely strong disorder. This study is motivated by 



two reasons. First, the Efetov-Tschersich linear law 
for the conductivity is by no means identical to the ex- 
perimentally observed linear law Q for resistivity. It 
is important to understand, if the deviations from the 
law @ due to the disorder can explain this discrepancy 
between theory and experiment. 

Secondly, one can expect inhomogeneous fluctuations 
of gij to be indeed very strong. The intergrain conduc- 
tance gij exponentially depends on the thickness dij of 
an insulating layer between the grains 

9ij = a^ hi3 1 h ij = -2/e(dy -(d)), (5) 

where g = exp(lng). The decrement of the electronic 
wave- function within the insulating layer k ~ 1 — 5 A , 
so that already fluctuations of dij with quite moderate 
standard deviation of order of few Angstroms may cause 
gigantic fluctuations of gij. Thus, the condition of strong 
fluctuations 

a = » i (6) 

can be fulfilled very easily. 

Critical networks. Under the condition @ neither 
the perturbation theory, nor the effective medium ap- 
proximation, used, respectively, in the case of weak and 
moderately strong disorder in [8|, can provide a reliable 
tool for finding the conductivity of the system. One can 
use instead the theory of highly inhomogeneous media 
(see, e.g., 0,0])- According to this theory, the conduc- 
tivity of highly inhomogeneous network of conductances 
is determined by the so called critical subnetwork, con- 
structed according to the following recipe (we will call it 
the "g-procedure" ) : 

• The critical conductance g c = ge ha is found from 
the condition 

/>oo 

/ P(h)dh= Pc , (7) 
Jh c 

where P(h) is the distribution function of h (we as- 
sume that there is no correlation between conduc- 
tances of different contacts), p c is the percolation 
threshold for the "bonds problem" on a lattice of 
granules (regular, or irregular - depending on the 
geometry of the array). 

• Certain number q is chosen, such that 1 ^ 9 « A. 

• All the intergrain contacts with hij < h c — q arc 
substituted by infinite resistances (disconnected). 

• The contacts with hij > h c + q are substituted by 
ideal conductances (abridged). Thus, the entire set 
of granules falls apart into a number of finite clus- 
ters, all members of the same cluster being ideally 
connected. These clusters serve as vertices n in the 
effective network. 
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FIG. 1: A fragment of a critical network C) 1 J 1 7 , corresponding 
to a particular random conductance network with A = 8.66 
on a square lattice 140 x 140. The vertices (clusters) are shown 
as black circles, a thickness of a bond between two vertices 
n, n' is proportional to the connection number N nn i . 

• The "critical conductances" with \h — h c \ < q can 
cither connect granules within the same cluster n 
(then they are shunted by ideal conductors and 
can be discarded), or establish links between dif- 
ferent clusters n and n' . Thus, we arrive at an 
effective network of vertices n, some of them be- 
ing connected by critical conductances. There may 
be more than one critical conductance, connecting 
the same pair of vertices. By construction (be- 
cause h c + q > he), besides finite fragments, there 
must be also an infinite connected subnetwork in 
this effective network, spanning through the entire 
system. As usual, one can cut off all "dangling 
ends" of this infinite subnetwork and thus extract 
the backbone, which is the substructure, relevant 
for the dc conductivity (for the consistent defini- 
tion of the backbone and the dangling ends, see, 
e.g., [13 )■ It i s this backbone, that we will call the 
d-dimensional critical network Cq of width q. A 
fragment of a typical computer-simulated critical 
network is shown in Fig^ 

Finally we are left with a random multiply connected 
network, in which every conductance is an independent 
random variable. The geometry of this network is highly 
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FIG. 2: The probability distributions for coordination num- 
bers z and connection numbers N for the same critical net- 
work, as in FiglTI 



nontrivial. Distribution functions for coordination num- 
bers z n , and for the "connection numbers" N nn r (num- 
bers of critical conductances, connecting vertices n and 
n') are shown in FigEI It should be noted that z n at 
different vertices, as well as N nn / at different links, are 
correlated, because vertices, corresponding to larger clus- 
ters of abridged grains, have a tendency to both larger 
coordination and larger connection numbers. 

The statistical geometry of the critical network, no 
matter how complicated, is universal: it depends neither 
on details of the initial lattice of grains, nor on the distri- 
bution function of individual conductances. Topological 
properties of critical networks for all systems of a given 
spatial dimension d are statistically equivalent. Different 
critical networks within the same universality class differ 
from each other only in the spatial scale - the correlation 
radius (see (lif): 



= a[P(h c )Y 



iA- 



(8) 



where a is the typical intergrain distance, v& is the criti- 
cal exponent of the correlation radius for the percolation 
problem in d dimensions (yi = 4/3 and w 0.875, see 
[ljj). The conductivity of the critical network is 



2ttH q y 
.(<i) 



2irha d - 2 



(9) 



where A q is a numerical prefactor. 

Since the critical network C q explicitly depends on 
the value of the parameter q, so does the corresponding 
conductivity o~ q . To get physically relevant o~q—>oo one 
has, in principle, to extrapolate the result to q 3> 1. For- 
tunately, already for q > 1 the value of q affects only the 

numerical preexponential factor Aq , which rapidly con- 
verges to the true value A^ as q grows (the corrections 



are exponentially small at q 3> 1). The conductances 
with | hij — h c \ 3> 1 do not play any considerable role in 
the critical subnetwork, therefore for a qualitative dis- 
cussion it is sufficient to choose some q ~ 1 and not to 
bother much about the effects of finite q. These weak 
and purely quantitative effects will be discussed below, 
when we will compare the numerical results, obtained for 
different values of q. 

Temperature dependence of the critical net- 
work's conductivity. In our problem the conductiv- 
ity is also determined by the critical network, but the 
conductances in this network evolve with "RG time" t 
according to the RG equations 



dt 



(10) 



where index k = 1, . . . , N nn > distinguishes different con- 
ductances, connecting the same pair nn' of vertices; R nn ' 
is the resistance of the critical network between vertices 
n and n' . It is very important, that these resistances are 
governed solely by the critical conductances, so that the 
system (fTUf) is closed. 

Let us introduce normalized variables 



t 



A (d h 



AW, 



d - 2 2irh 

(ii) 



then 



• The critical network is reduced to a completely uni- 
versal one, with the spatial scale equal to unity. 

• The equations © take the universal form 



dy. 



(fe) 



^^({e""}). (12) 



Initial conditions y^/(t = 0) for the RG equations 
l|12|l independently take random values with a uni- 
versal distribution function 



fl(g-M) 

2q 



(13) 



Thus we have arrived at a completely universal math- 
ematical problem: One has to solve a set of universal 
equations (with universal initial conditions) for variables 
Vnn' ' defined on a random graph with universal statisti- 
cal properties. As a result, the normalized conductivity 
a should be a universal function of the normalized time 
t: 



(14) 



the shape of this function depends only on the space di- 
mensionality d and on the parameter q. At q > 1 the q 
dependence is very weak and can be neglected, therefore 
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we omit the index q in what follows. Unfortunately, we 
were not able to find this universal function analytically, 
so we had to do that numerically (see below). 

Coming back to the initial physical variables, we get 
the global conductivity of the system in the form: 



p(d.) 



e 2 /2irh 



Tr 

[^)] d " 2 a(T«) m T 



In 



(15) 



where T R = g c Ec- The only nonuniversal (i.e., depend- 
ing on the details of the conductances' distribution func- 
tion) ingredient in this formula is £W ~ A" d . Note also, 
that there is absolutely no reminiscence of the nonuni- 
versal short-scale geometric structure of the array. 

Lateral conductance of granular films. One 
should distinguish between thin films, whose thickness 
SL « £ (3) , and thick films with SL > The ar- 

guments of paper (see also book 03 )j designed for 
the description of hopping conductivity in films, can be 
modified for our case. The true 3-dimensional behavior 
takes place only in thick films, where the dimensionless 
conductance of a square sheet of the film is Go — cr^SL, 
and 



G D (T) 



G a (T R ) 
For thin films 

Go(T) 



= f& 



SL 



G a (T R ) 



jr(2) 



HTr/T) 
G D (T R ) 



G a (T R ) 

and Gu exponentially depends on SL: 

G a (T R )=AWg Q e h °( 5L \ 



(17) 



(18) 



where h c = h c (SL) is determined as a threshold of lateral 
percolation in a slab of thickness SL cut out of a three- 
dimensional network of all conductances with hij < h c . 
The threshold condition has a form 



^ 3 \h c ) = cSL, 



(19) 



where the three-dimensional correlation length £( 3 ) (h c ) ~ 
a[(h c — /i c )/A] _l/3 , and c ~ 1. As a result, in the range 
a < SL < £ (3) one finds 

h { c 3) - h c {SL) = DA{a/5L) 1/u3 , 



(20) 



where D is some nonuniversal numerical constant (see 
[lif). Note, that this formula interpolates between 

(2) 

in the case of a bulk, thick film, and hi ' in the case of a 
one-monolayer film. 

~(2) 

The quasi two-dimensional critical network C q , aris- 
ing in the film as a result of the g-procedure, centered at 
h c (SL), is characterized by the effective two-dimensional 

(2) 

correlation length . To estimate this length we write 



it in a form £ 



(2) 
eff 



the scaling arguments. This expression, valid in the 
range a <C SL <C should match with the true two- 



dimensional behavior £ 



(2) 
eff 



^(2) _ A" 2 at SL ~ a, and 



with the true three-dimensional one Q s = ^ 3 - ) ~ A 1 ' 3 
at SL ~ £( 3 ) . These two matching conditions fix the ex- 
ponents ai,aa, so that we find a\ = 1 — (^2/^3) and 
«2 = ^2- As a result, can be conveniently written in 
a form 



t(2) 
?cff 



SL 



(f^/SL) »5L. 



(21) 



a(8L/a) ai A" 2 , which follows from 



Thus, our quasi- two-dimensional critical network is in- 
deed effectively strictly two-dimensional, hence it belongs 
to the corresponding d — 2 universality class, and the 
universal function is F^> for all films with SL <C S,^ (cf- 
formula j!7|)). not only for those, consisting of just one 
monolayer of grains. 

Numerical results. In this Letter we present the nu- 
merical results only for the two-dimensional case. For 

(2) ~ 

determination of the function Fq (t) we have gener- 
ated random L x L square arrays of conductances with 
distribution P(h) = 0(\/3A - \h\)/2y/3A for different 

20 < L < 140 and < A < 8.66. For the samples 

(2) 

(16) with large A we constructed critical networks C q with 
0.5 < q < 2, and numerically solved the system of equa- 
tions I jl 2(1 and (@J on these networks, which are very much 
smaller arrays, than the initial ones. Note, that only due 
to this truncation of the array, as well as to the fact, 
that the conductances of the critical network are all of 
the same order of magnitude, the problem becomes fea- 
sible. The numerical solution of the full set of Kirchhof 
equations for a reasonably large array of strongly differ- 
ent conductances would be practically impossible. Thus, 
the concept of the critical network turns out to be cru- 
cially helpful not only for understanding of basic physics, 
but also for the numerics. 

The main result - the dependence of the conductivity 
on \n(T R /T) for different values of A is shown in FigOU 
One can see that this dependence becomes universal al- 
ready at A ~ 3 — 4: the difference between curves with 
A = 3.47 and A = 8.66 is negligible, so we can with 
reasonably high accuracy identify the last one with the 
universal function F^(t). In the range < t < 0.5 it is 
approximately linear with a slope, corresponding to an 
effective coordination number z c g w 3.2, which is consis- 
tent with the statistics of coordination, shown in Fig|2] 
Note, that z c g is universal and does not depend on the 
coordination number z of the initial lattice. At t > 0.5 
the slope starts to decrease rapidly, and F^(t) goes to 
zero at t c w 1.2. Unfortunately, the problem of finding 
the critical exponent, which characterizes the behavior of 
F^(t) at t — > t c , seems to be numerically not feasible. 

We also have checked the accuracy of the method from 
the point of view of finite size and finite q corrections. 
For an array 40 x 40 with A = 1.734 we were able to 
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FIG. 3: Dependence of the normalized conductance per 
square for a thin film on the parameter t = \a(Tn/T) for dif- 
ferent standard deviations A of the logarithms of intergrain 
conductances. M is the number of samples in the averag- 
ing ensemble. In all cases q = 1.7 and L > £ (2) (A). (a) 

- A = 0.14, L = 50, M = 150; (b) - A = 0.58, L = 50, 
M = 150; (c) - A = 1.16, L = 50, M = 300; (d) - A = 2.31, 
L = 80, M = 300; (e) - A = 3.47, L = 100, M = 400; (f) 

- A = 8.66, L = 140, M = 400. in all cases q = 1.7 and 
L » £ (2) (A). The thick line - universal function F {2) (t). 



perform the straightforward calculations (involving the 
entire array) and compare their results to those, obtained 
on truncated C^-networks. The results of this compari- 
son, shown in Fig^J demonstrate a very good accuracy 
already at moderate values of q = 1.7. This accuracy 
is reached at the network, which includes only about a 
half of all conductances of the array. Note, that although 
the condition q <C A, required for general distributions 
P{h), is obviously not satisfied in this particular simu- 
lation, for the rectangular shape of P(h), used here, the 
strong inequality can be substituted by a milder condi- 
tion q < V3A. 

For a sample 140 x 140 with A = 8.66 (the same realiza- 
tion, as used in Figs lll2J) we performed the calculations 
with different values of q. The results, shown in Fig0 
demonstrate fast convergence of the procedure: the dif- 

(2) 

ference between Fq for q = 1.7 and q = 1.9 is already 
very small, so that using q — 1.7 procedure for evaluating 
the universal function is a very good approximation. 

From general arguments it is clear that to get self- 
averaging results, one should use arrays with size L > 
£( 2 ). In FigEJwe have plotted functions F^ 2 \t), found for 
arrays with the same A = 8.66, but with different sizes 
L. For these arrays one expects £( 2 ) ~ C[P{h c )]~ U2 — 
C ■ 30 4 / 3 = 93C (with unknown numerical coefficient C). 
We see that for L > 100 the curves become almost iden- 
tical, so we conclude that C ~ 0.5 — 1 and square arrays 
with L > [P(/i c )] _!y2 can empirically be considered as 




t/a(0) 



FIG. 4: Comparison of the solutions of the full set of equation 
for the 40 x 40 sample with A = 1.734 (solid line) and the set, 
truncated according to the q-procedure with q = 1.7 (broken 
line) . 




FIG. 5: Comparison of the solutions, using (/-procedure with 
different q. One can see that the difference between q = 1.7 
and q = 1.9 results is already negligible. 



large enough. 

In principle, the proper thing to do would be an accu- 
rate extrapolation of fmite-q and fmite-L results to q = oo 
and L — oo. Unfortunately, the complexity of the nu- 
merical solution at large A, q, and L did not allow us to 
reach the accuracy, necessary to such an extrapolation. 
The larger q is chosen for g-procedure, the less stable and 
precise the calculation of Rij becomes. This instability 
limits q by about q = 2. Another precision limitation is 
rapid increase of computational complexity with increase 
of L/t;( 2 \ The mentioned difficulties limit the accuracy 
of each curve by about 5%, which is comparable with the 
separation of neighboring curves on our plots, making 
the extrapolation procedure senseless. However, since a 
rapid convergence of the results is observed, we believe 
that the universal curve F^ 2 \t) lies within the same 5% 
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FIG. 6: Comparison of solutions for different lattice sizes L. 
All curves are obtained for largest standard deviation A = 
8.66 by the g-procedure with q = 1.7, and averaged over M = 
400 realizations. 

error interval near the curve with A = 8.66, L = 140. 

Conclusions. Thus, the dependence of conductivity a 
on In T for strongly disordered arrays of grains is nonlin- 
ear and universal, contrary to the case of a regular array, 
for which this dependence is linear, with a nonuniversal 
slope. The sign of the curvature of this nonlinear func- 
tion is such, that the effects of strong disorder seem to 
improve the agreement between the theory and the ex- 
perimental law Q). Unfortunately, the available experi- 
mental data is still not sufficient for a direct quantitative 
comparison of theory and experiment. In particular, it is 
hard to make reliable estimates for Tr and Gd(Tr). It 
is also worth reminding that the one- loop RG, which is 
the basis of the present theory, is valid at temperatures 
well above the metal-insulator transition. 

The universality of the temperature dependence of the 
conductivity is the direct consequence of the universality 
of the critical network, which carries the current in a 
highly inhomogeneous system. This critical network in a 
granular film appears to be effectively two-dimensional, if 



the film's thickness 8L is less than the three-dimensional 
correlation length £^ 3 ) ~ aA U3 , the latter being large at 
high level of disorder. 
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